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ON THE INTEGRABIUTY CASES OF THE EQUATION 


OF MOTION FOR A SATELLITE IN AN AXIALLY 
SYMMETRIC GRAVITATIONAL FIELD 


Abohghassem Ghaffari 


ABSTRACT 

The projection of an axially symmetric satellite’s orbit on a plane perpen- 
dicular to the rotation axis (z = const.) is given by the second-order differential 
equation 


-IL = *, - y' . (i) 

1 + y 

where the prime denotes the derivative with respect to x and \p (x, y) is a known 
function. Two integrability cases of (1) have been investigated and it has been 
shown that for these two cases the integration of (1) can be carried out either by 
quadratures or reduced to a first-order differential equation. Analytical and 
physical properties of (1) are expressed, and it is shown that equation (1) can be 
derived from the classical plane eikonal equation of geometric optics. 


r -rni\'G PA^c 


vt 



SUMMARY 


This paper discusses the integrability cases of a special equation of motion 
recently presented by H. Knothe* in his article "Satellites and Riemaniann 
Geometry". The importance of this equation (equation 1) is its close connection 
with the plane eikonal equation of geometric optics. The equation under dis- 
cussion is given by 


-Z-— = - y' . (i) 

1 + y' a ? 


y(x) represents the projection of the satellite's orbit on a plane z = const., 
prime denotes the derivative with respect to x, and 0 (x,y) = log 4' (x,y) , where 
\p (x,y) is being interpreted* as the speed with which the projection of satellite 
moves in the two-dimensional space (x ,y) . 

Knothe* has shown that for a given potential 0(Uj, u 2 ) in an axially symmetric 
three-dimensional Riemaniann space (Uj, u 2 , U 3 ) there exists a function 
v//(u , u 2 ) by means of which the equations of motion for all point- masses with 
equal total energy and equal angular momentum can be reduced to a single 
ordinary second-order differential equation. Considering the special case of 
the motion of a satellite around a rotationally symmetric body in three- 
dimensional Euclidean space (x,y,z), Knothe* has shown that the projection of 
a satellite's orbit about a body of revolution on a plane orthogonal to the axis of 
rotation (z-axis) is a geodesic curve given by equation (1) . Equation (1) can also 
be considered as the autonomous equation of motion of a point-mass moving in 
a field of force with a force function ^ (x,y) = e^ ( x ’ y * • 

In this paper two integrability cases of equation (1) are considered and it is 
shown that for these two special cases the integration of ( 1 ) has been reduced to 
quadratures and to an ordinary first-order differential equation. Analytical and 
physical interpretations of this equation as well as its connection with the classical 
plane eikonal equation of Geometric Optics are also expressed. 


•Celestial Mechanics 1 (1969), PP- 36-45; D. Reidel Publishing Company, Dordrecht-Holland. 


ON THE INTEGRABIUTY CASES OF THE EQUATION 


OF MOTION FOR A SATELLITE IN AN AXIALLY 
SYMMETRIC GRAVITATIONAL FIELD 


Abohghassexn Ghaffari 


1. Introduction 

This paper has been devoted to the investigations of the two integrability 
cases of the ordinary second-order differential equation 


-2 — = </■„- y' ^ , (i) 

1 + y' 2 y 

which defines the projection (on a plane z = const.) of the motion of a satellite 
around an axially symmetric body in the three-dimensional Euclidean space 
(x, y, z ). Prime denotes the derivative with respect to x. 

Recently H. Knothe (ref.) showed that for a given potential 0(u lf u 2 ) in an 
axially symmetric three-dimensional Riemannian space (u j , u 2 , u 3 ) there exists 
a function 0 (u t , u 2 ) by means of which the equations of motion, for all point-masses 
with equal total energy and equal angular momentum, can be reduced to a single 
ordinary second-order differential equation. 

Considering the special case of the motion of a satellite around a rotationally 
symmetric body in three-dimensional Euclidean space (x, y, z) , Knothe (ref. p, 45 , 
formula 65) has showed that the projection of satellite's orbit on a plane z = const. 
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is a geodetic curve given by the equation (1) , where 


^ - log (x, y). 


( 2 ) 


The function 0 (x, y) has been interpreted by H. Knothe (loc. cit.) as the 
speed with which the point-mass moves in a two-dimensional Euclidean space 
(x, y), i.e., 

* = 37 <3> 

where a is the arc length of the orbit's projection on the plane z = const. 
Equation (1) is a special case of the general equation 


( z l + z l + 1 )y" ~ z 


z y 

X yy J 


•3 


+ (2 z z -z z 
v X xy y yy 


)y ,2 + (z 2 -2z z )y'-z z 

/ j \ x yxy 7 * 7 yxx 


( 4 ) 


giving the geodetic curves on a given surface 


z = z (x, y) . 


The second-order equations (1), (4) are of the third degree in y' and there are 
no general rules for solving nonlinear differential equations, and skill and ingen- 
uity are essential. The purpose of this paper is to show that there are (besides 
the trivial case in which the variables in I p can be separated), at least, two cases 
in which the integration of equation (1) can be carried out either by quadratures 


Z 


or reduced to an ordinary first-order differential equation. These cases are 
deduced from the form and structure of the function 0 expressed in both rec- 
tangular and polar coordinates. 

2. Integrability Cases 

a. Equation (1) can be written 

dv # i — 

- d arct y = 0 dx - 0 dy . (5) 

. 12 y * 

1 + y 1 

The integration of (5) would be easy if the right-hand side of (5) is also an exact 
total differential of a function U (x, y) such that 


0 dx - 0 dy = dU = U dx + U dy • 

/ * * y 


(6) 


The integrability conditions of (6) leads to the Cauchy-Riemann differential 
equations 


r 



provided 0 e C 2 and UeC 2 in a domain D of (x, y) of the space R 2 . 

The infinite number of functions U (x, y) are given by the two following 
quadratures 


I 


U (x, y) = I 0 y (x, y) dx 




0 X (x 0 , y) dy + a 


( 8 ) 


where x 0 and y 0 are arbitrary numerical values of x and y respectively, and a 
is an arbitrary constant. The particular function U corresponding to a = 0 is 

y 

(x 0 , y) dy, (9) 

0 

and the integration of equation (5) gives 


U (x, y) = 


/ 

•'X , 


0 y (x, y) dx 


-I 


arctg y' = U (x, y) + C, y* = tan (U + C) = tan oj 


or 


U + C = w 


(10) 


where the angle has been defined by H. Knothe (ref.) such that the tangent vector 
of the trajectory x(t), y(t) has the components cos co and sin to respectively, and 
C is an integration constant. 

Analytic and physical inteipretations . 

The Cauchy-Riemann conditions (7) state that there exists an analytic 
function W of the complex variable z = x + iy such that 


W (x + iy) = U (x, y) + i 0 (x, y) , (11) 

where the two harmonic functions U and 0 satisfy Laplace equations 


AU = 0, A 0 = 0. 


( 12 ) 


From the Cauchy-Reimann conditions (7) one can also deduce that 


U 2 + IJ 2 = 3 2 + d, 2 

x T y r x r y 


or the expression 


is an invariant (in terms of x and y) on either surfaces 


( 13 ) 


(14) 


r z s "J (x, y) 

< _ 

* - (x. y; • 


(15) 


There is a close relation between equation (1) and the classical plane eikonal 
equation of geometric optics 


(Vz) 2 = z 2 + z 2 = n 2 (x, y) (16) 

where n is a function of x and y. In fact, equation (1) can be derived, for a special 
form of n (x, y), from equation (16). The form of equation (16) suggests to 
write 

n (x, y) = e* , 




(17) 


and setting 


z e* cos u) 


z _ sin w, 


(18) 


one finds that u> and o> satisfy the first-order partial differential equation 

x y 


cos co • to r s in co • co = 0 cos co - 0 s in 

x y ~y r x 


(19) 


with the corresponding characteristic differential equations 


, dy da> 

dX = 1 SB" 

t g W 0 y - 1 ,* OJ 


( 20 ) 


which give 



dw , -7 -7 

tg a), and " = 'Py ~ ’/'x 


tg 


( 21 ) 


Taking the derivative of the first equation of (21) with respect to x and 
considering the second equation of (21) we get exactly the equation (1). 

The case n = const, corresponds to the developable conoid surfaces whereas 
n as a function of x and y corresponds to the plane wave fronts of light in an 
inhomogeneous medium with a variable index of refraction n (x, y). 


If n is a function of z alone, equation (16) can be integrated by classical 
methods and represents the well-known parallel (moulding) surfaces. 

Equation (1) can also be considered as the plane autonomous equation of 
motion of a point-mass moving in a field of force with a force function 

0( x, y ) = e ^ x>y) . 

b. If 0(x,y ) is an homogeneous function of degree zero, that is to say if 

'P (x, y) = 0 ^1, L.) - f = f (u) (22) 

where 


u 



(23) 


then the equation (1) will be reduced to an ordinary first-order differential 
equation. The new variables are x and u. 

In fact, equation (1) becomes 


= f' (u) (1 + u y' ) , 

1 + y' 2 


and the substitution (23) gives 


y' =u'x+u = t + u 


^u'x + 2u'= u' l^- + l) 


\ ax ) 


and equation (24) becomes 


t ^ + I'j = f ' (u) [1 + u (t + u)] [1 + (t + u) 2 ] 


which is of the first order in u and t. 


(24) 


(25) 


(26) 
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If, in polar coordinates, 0 is a function of the argument 

9 = arctg^ 
x 

alone, the substitution (23) leads to the same equation (26). 

It may exist, besides the two above mentioned cases, other ca^es in which 
the integration of (1) can be carried out by quadratures or reduced to an ordinary 
first-order differential equation. 
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